Invariant correlational entropy as a signature of quantum phase transitions in nuclei 
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We study phase transformations in finite nuclei as a function of interaction parameters. The 
signature of a transition is given by invariant correlational entropy that reflects the sensitivity of an 
individual many-body state to changes of external parameters; peaks in this quantity indicate the 
critical regions. This approach is able to reveal the pairing phase transition, identify the isovector 
and isoscalar pairing regions and determine the role of other interactions. We show the examples of 
the phase diagram in the parameter space. 
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Recently, an appreciable effort has been applied to un- 
derstand and classify quantum phase transitions "7, '5, "s*! . 
Although the concept of a phase transition is strictly ap- 
plicable in the thermodynamic limit only, there are nu- 
merous examples of structural changes in mesoscopic sys- 
tems, ranging from molecular clusters and semiconduc- 
tors to atomic nuclei or quark systems, under variation of 
control parameters. Proper understanding of such tran- 
sitions is crucial for areas from cosmology and formation 
of the universe to quantum computing and decoherence. 
The typical feature of phase transitions in small systems 
is the absence of discontinuities in the observables, and, 
therefore, difficulty for identification and classifications of 
such transitions. The counterparts of phase changes in 
small systems involve restructuring, critical sensitivity to 
parameters, and chaotic large-scale fluctuations. In this 
work we study mesoscopic phase transitions in atomic 
nuclei using shell model interactions, which are known to 
reproduce the low-lying states in selected nuclei with a 
remarkable quality. The instrument we suggest for such 
studies can be similarly applied to other finite quantum 
systems. 

We identify the presence of a phase transformation as 
an enhancement in the invariant correlational entropy 
(ICE) which was introduced in 4] . ICE provides a mea- 
sure of sensitivity of a given state in the many-body sys- 
tem to variations of external parameters. From earlier 
studies of bosonic models the increase of ICE is known 
to be associated with critical points 0. We do not con- 
sider here thermal phase transitions for a system in a heat 
bath although it is possible, see 3 and 0, to introduce 
an effective temperature as a measure of fluctuations in 
the dynamical response to the change of parameters. 

Following 01, we assume a Hamiltonian H{X) that de- 
pends on a parameter A, so that, in an arbitrary basis 
|fc), any eigenstate |a(A)) of H{X) can be decomposed as 
|a(^)> = Efc CkW\k) ■ The ICE is then defined as 

5"(A) = -Tr{p"ln(p")}, (1) 

where is the density matrix of the state \a) in the 
basis k averaged over a small region A G [A, A + 5], 



The discussion of various entropy-like quantities of in- 
dividual wave functions can be found in Ref. ^. ICE 
is basis-independent von Neumann entropy that reflects 
the correlations between the wave function components 
along the evolution path determined by the parameter A. 
It should not be confused with basis-dependent Shannon 
information entropy 



(3) 



that was extensively used for studying the degree of com- 
plexity of the state \a) with respect to the given basis \k) 
0, 01 without taking into account correlations between 
the amplitudes C^. 

For a pure state, i.e. without averaging, the density 
matrix has a single non-zero eigenvalue equal to one be- 
cause of the normalization Trp = 1. This results in 
S"" = 0. After averaging, we deal with a mixed state, and 
the eigenvalues of p deviate from the trivial limit. As 
shown in higher orders of perturbation theory bring 
in new non-zero eigenvalues. Averaging over only two 
discrete points A and A' leads to a factorized matrix 
with two non-zero eigenvalues 



4 = (l±|(a(A)KA'))|)/2. 



(4) 



Now ICE can change from (in the absence of any evo- 
lution of the state) to the maximum value In 2 reached 
for the orthogonal states in the representative points A 
and A'. Thus, ICE in a basis-independent way shows 
how "quickly" the eigenstate a reorients due to sudden 
change of interaction from H{X) to H{X'). This makes 
ICE an ideal tool for studying phase transitions. 

We apply the ICE tool to spherical shell-model sys- 
tems, where the general Hamiltonian includes indepen- 
dent particle orbitals and two-body residual interactions. 
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Here pair creation and annihilation operators couple cor- 
responding single-particle operators to an appropriate 
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angular momentum, i.e. pf^^) = (H-(5i2)~"^''^[ai x a2]LA. 
The usual pairing interaction is given by the L = 
term. We suppress the isospin degree of freedom, but 
it is straightforward to include it in direct analogy to 
angular momentum. 

We start with a simple model, where the pairing in- 
teraction is known to result in a phase transition in the 
BCS approach. N identical nucleons occupy two single- 
particle orbitals ei = — e and 62 = e of equal spheri- 
cal degeneracy ili — 0.2 — fi. They interact via off- 
diagonal pair transfers determined by a parameter A as 



^Jii;22) ^ y(^^-^'-^i ^ -2X/n. Let the system be half- 
occupied, N = il, which sets the chemical potential 
/Lt = 0. It is known from the BCS theory that the pairing 
interaction causes the normal to superconducting phase 
transition, and, in the asymptotic limit of a large system, 
any weak but attractive pairing is sufhcient for creating 
the Cooper instability. For a finite system the BCS solu- 
tion is approximate. Because of the finite single-particle 
level spacing at the Fermi surface, the BCS phase tran- 
sition takes place at the critical strength value. 



^(22:11) 
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The exact solution of the problem [3 determines the 
ground state and the excited states of two types, "pair- 
vibrational" states with seniority s = and redistribu- 
tion of the pairs, and states with broken pairs and s ^ 0. 
For each state one can calculate ICE varying A. The re- 
sults for the ground and two lowest pair-vibration states 
calculated with the averaging interval 6 — 0.01 are shown 
in Fig. n] The clear ICE peak in the ground state in- 
dicates the presence of the phase transition. For e = 1 
and n ~ N — 16 the BCS theory predicts a phase transi- 
tion at Ac = 0.25. Corrections to the BCS vacuum [1,0 
make Ac slightly larger, and ICE reaches its maximum 
at Ac = 0.3. In the excited states the phase transition 
is quenched. The first excited state is collective and ex- 
hibits a broadened peak of ICE. For the second excited 
state the peak disappears although there is still a maxi- 
mum of S{X) at A = indicating an enhanced sensitivity 
to switching on pairing. 

It is also possible to recognize the signature of the 
phase transition in the exact energy spectrum. In Fig. 
121 the excitation energy of the first pair vibration state 
is shown. The softening of this mode at the same point 
Ac = 0.3 is characteristic for the phase transition and sup- 
ports the level crossing interpretation. With the dashed 
line we show the same state in the random phase approx- 
imation, RPA, built on the BCS condensate The in- 
stability of the BCS-I-RPA approach for a soft mode is 
an artifact of the invalid approximation, see also 01 ■ 

As a realistic example we consider the case of the 
TV = 28 subshell closure in ^^Ca. It has been empha- 
sized in Ref. ^3 that although the BCS approach pre- 
dicts no paired ground state in this nucleus, the exact 
solution indicates presence of pairing correlations. The 
interplay of weak pairing and self-consistent monopole 
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FIG. 1: Invariant correlational entropy of the ground and two 
lowest pair-vibrational excited states calculated for the two- 
level pairing model as a function of pairing strength A and 
with averaging interval S — 0.01. 
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FIG. 2: Energy of the first pair-vibration level in the 
BCS-I-RPA (dashed line with instability at Ac) and exact so- 
lution (solid line) that shows softening at the shifted point. 



renormalization of single-particle energies leads to corre- 
lation energy of nearly 2 MeV. The specifics of pairing 
correlations and self-consistent mean field in TV = 28 nu- 
clei have recently been a subject of experimental [T]| and 
theoretical 0,0] interest. Here we go beyond the mean 
field and solve exactly the many-body Hamiltonian with 
the effective interaction |0 in the plane of parameters 
Xp and Xnp that scale all pairing (T = 1, L = 0) and 
remaining non-pairing matrix elements, respectively. 

The results for the ground state of '*^Ca are given in 
Figs. 13 and 2] The peaks in ICE in the upper panel in 
Fig. 13 unambiguously identify the pairing phase transi- 
tion. The drop in the occupancy of the lowest /y /2 orbital 
shown in lower panel is also consistent with this behav- 
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FIG. 3: ICE (computed with averaging interval S = 0.1), 
upper panel, and occupancy of the lowest f-r/2 orbital, lower 
panel, as a function of pairing strength for different scales of 
non-pairing interactions, Xnp =0.2, 1, 2, 3, and 4. 



ior indicating the disappearance of subshell closure and a 
transition into a paired state. Five curves in both panels 
of Fig. O allow one to track the evolution of the phase 
transition as a function of strength A„p of non-pairing 
matrix elements. Panel (a) in Fig. 01 shows ICE with 
a contour plot in the Ap - A„p plane. Higher values are 
indicated with shaded areas. ICE is enhanced along the 
diagonal outlining, as in a phase diagram, the bound- 
ary between the normal Fermi state (lower right) and 
superconducting paired state (upper left). The panel (c) 
confirms this identification of phases. Here the effective 
number of pairs Mt=i , computed as an expectation value 
of the operator 
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(pill) y 
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in the ground state, is shown. This quantity reaches its 
maximum value A/i — {1/A){N — s){^l — N — s + 2), where 
f2 is the total capacity of the valence space, and seniority 
s gives the number of unpaired particles, in the degener- 
ate limit, i.e. when the pairing becomes so strong that 
differences in single-particle energies can be ignored. The 
maximum number of pairs is expected for a fully paired 
state with 5 = 0, which for our model, with ~ 20 
and N — 8, leads to A/i — 28. The non-pairing interac- 
tions create a random background, mixing seniorities and 
adding some statistical number of pairs to any state. This 
background can be estimated by averaging operator A/i 
over aU many-body states, J/^ = [N{N - l)]/[2(fi - 1)], 
which gives A/7 « 1.5 in our model. Indeed, the major 
part of the valley in Fig. 2{c), shown as a white area. 



is levelled at around this value. In the phase transition 
region, along with the enhancement of ICE, the number 
of correlated pairs quickly rises, finally becoming close to 
28. The full saturation at A/i — 28 is not expected to be 
reached, since it is only possible for the constant pairing 
interaction. 

The enhancement of ICE not only allows one to lo- 
calize the phase transition but also quantifies the sharp- 
ness of the transformation and the size of the critical 
region. In macroscopic BCS theory it is usually assumed 
that the non-pairing interactions renormalize quasiparti- 
cles but essentially do not participate in the phase tran- 
sition. The shell model analysis demonstrated that 
this is not the case, at least in a finite system. As seen 
from Figs. 13 and 0] the scaling factor A„p plays a signif- 
icant role in the overall picture. At A„p = 0, the pairing 
phase transition is quite sharp and takes place slightly 
below Ap = 1.0. Indeed, the comparison of earher studies 
[Tol ITsI shows that pairing would be stable and treatable 
with the BCS if no other interactions had been present 
in ''^Ca. Presence of other parts of residual interaction 
softens and widens the transitional region, putting the 
physical state A„p = Ap = 1 in the normal domain below 
the phase transition peak, although within the region of 
enhanced pair fluctuations. It was also shown earlier 
that other interactions destroy the purity of the seniority 
classes (in the first order the ground state acquires the 
contribution of s 4) moving the dynamics to many- 
body quantum chaos (tI, Q . 

The non-pairing interactions shape the properties of 
the mean field and redefine effective single-particle en- 
ergies; on the other hand, they induce chaotic residual 
dynamics. In Fig. E^b), the occupancy of the fj/2 or- 
bital is mapped in the same parameter space. The lowest 
fj/2 level in ^^Ca is separated from other single-particle 
levels by almost 2 MeV. At weak residual interactions 
the ground state is the slightly perturbed Fermi sea, 
with nearly all A^^ = 8 valence neutrons occupying the 
fj/2 level. This region is seen as the dark area in the 
lower left corner in Fig. EJb). As Ap grows, the pairing 
phase transition breaks the magicity, and the depopu- 
lation of the orbital is consistent with the line of 
pairing phase transition, see also Fig. |3| Along the A„p 
axis, the original state with the fully occupied lowest or- 
bital also disappears around A„p = 4, similarly to the 
change of magic numbers far from stability. The bare 
single-particle energies correspond well to the mean field 
defined with respect to the ^"Ca core. The increase of 
residual interactions renormalizes the mean field and at 
some point makes the original single-particle basis incon- 
sistent with the potential. Strong residual interactions 
make the dynamics chaotic that can be interpreted as 
a non-zero single-particle temperature 0, . The occu- 
pation of the lowest orbital falls from A^7/2 = 8 to the 
"thermal" occupation A^7/2 — 8{N/U) — 3.2. When non- 
pairing interactions grow to around A„p = 5 and beyond, 
the pairing to normal phase transition again sharpens. 



Here the rise in ICE is much higher than at X„ 



0. 
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FIG. 4: Phase diagram of the shell model ground state wave function of Ca. The parameters are the pairing strength Ap and 
an overall strength Anp of remaining non-pairing matrix elements. Panel (a) shows ICE with the averaging interval S = 0.05 
in both scales. ICE is normalized by a factor 0.015. Panels (b) and (c) show the occupancy of the f-r/2 orbital and the average 
number of correlated pairs A/i, Eq. ||7J, respectively. To emphasize the consistency between these quantities and location of 
phase changes, the dotted contours in panels (b) and (c) outline the region where ICE is enhanced above the value of 0.5. 



In order to create a coherent paked state, pairing forces 
need to overcome not only the spread of single-particle 
energies but also random motion generated by strong 
residual interactions. 

From ICE in Fig. Efa) we learn that the evolution 
of the wave function driven by increase of non-pairing 
interactions is not followed by an enhancement of en- 
tropy. Since with a closed proton core a non-perturbative 
change in the mean field, such as onset of deformation, is 
not expected, the evolution of spherical orbitals, as well 
as chaotization of dynamics or rise in effective tempera- 
ture, are not associated with any phase transition. The 
role of A„p and the mean field in the pairing phase tran- 
sition in nevertheless eminent. Supporting earlier studies 
in [tI IT^IitI I. the sharpness of change and the transitional 
region between normal and paired states are influenced 
significantly by non-pairing interactions. 

As our last example, we consider the sd shell model for 
^^Mg. At N = Z, the competition of isovector, isoscalar, 
deformation and alpha-particle correlations makes the 
physical picture very complex. Einpirical information 
p^, studies of simple models [H El El El and 
analytic works (2^, supplemented by the direct shell 
model diagonalization [23, support the plausible exis- 
tence of different phases. However, it still remains un- 
clear whether these phases are pure enough to be iden- 
tified, or whether they are separated by distinct critical 
regions. Here the ICE provides a helpful tool. 

For this example we extend the list of parameters to in- 
clude Xt=i, At=o and A„p, which scale isovector (T = 1, 
L = 0), isoscalar {T — 0, L = 1), and all remaining two- 
body matrix elements in the Hamiltonian, respectively. 
The shell model is defined in this example with single- 
particle energies and interaction matrix elements from 



0.0025 



0.0020 



0.0015 



0.0010 



0.0005 



0.0000 




FIG. 5: The evolution of ICE: along the line of increasing 
isovector pairing At=i at fixed At=o = Anp ~ 1, solid line; 
as a function of At=o at fixed At=i ~ Anp = 1, dashed line. 
Interval S = 0.05 was used here. 



Ref. Eg. In Fig. El where ICE is shown as a function of 
Xt=i with Xt=o = A„p = 1 as a solid line and as a func- 
tion of At=Oj while At=i = A„p = 1, as a dashed line, in- 
dicates again a phase transition. The contour plot of ICE 
for this example is shown in Fig. We limit ourselves 
here to three different strengths of non-pairing interac- 
tions Xnp = 0.5, 1, and 2. This already gives a general 
idea of the ICE evolution. The valley in the vicinity of 
small At=o and At=i that corresponds to a normal state 
is clearly separated, and with increase of the strength of 
non-pairing interactions this region widens. However, at 
large A„p the sharpness of the phase transition is signifi- 
cantly reduced; note that the scale in the right panel in 
Fig. IHlis ten times smaller. This effect is analogous to 
the situation observed in ^®Ca, where, except for some 
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FIG. 6: ICE of the ground state in Mg computed on the plane At=o — At=i. The three panels show the evolution of the 
phase diagram depending on the overall scale of non-pairing matrix elements Anp- The averaging interval oi S — 0.05 was used, 
and the entropy is scaled (divided) by 0.015. 




FIG. 7: Effective average number of T = and T = 1 pairs, 
upper and lower panels, respectively, for the ground state of 
^*Mg at X„p = 1 on the At=o — At=i plane. The dashed line 
outlines the regions where ICE is enhanced above 0.45. 



special area, non-pairing interactions smeared the phase 
transition region. 

The only valley clearly separated by the mountains of 
ICE corresponds to the normal state. To identify phases 
of T = 1 and T = pairing, similarly to our previous 
example, we calculate the number of coherent pairs for 
A,ip = 1 case shown in Fig. [7| We consider isovector, 
T = 1, L = 0, and isoscalar, T — 0, L = 1, pairs 10. 



For both limiting cases of pure isovector and isoscalar 
pairing, A/r can be calculated using group algebra, see 
for example In the states of degenerate isovector 

pairing model 



1 



1 



-{N - s){n - iV - s + 6) + - 



(8) 



where t is isospin of unpaired particles, while T is total 
isospin of the system. Thus, in the state with maximum 
isovector pairing, s = 0, t = T = 0, one would expect 
TVi = N{n - N + 6)/8, which in our case of = 8 
and O = 24 leads to TVi = 22. Interchanging angular 
momentum and isospin and taking normalization of op- 
erators into account we obtain A/q — N{n — N + 6)/{2n), 
which results in maximum A/q ~ 3.7, whereas statisti- 
cally A/i « 1.8 and A/o ~ 0.3. This analysis along with 
Fig. [3 allow us to confirm the presence of isovector and 
isoscalar paired states in the upper left and lower right 
areas of the phase diagram. The isovector pairing phase 
appears as a shaded area in the lower panel of Fig. [T] and 
the shaded area of the upper panel indicates the isoscalar 
pairing state. Furthermore, the enhancement of ICE, the 
region outlined by a dashed line on both panels of Fig. 
[7| is consistent with the presence of the phase transition 
from normal to either of superconducting phases; this is 
the area of the most rapid rise in the number of coherent 
pairs A/o or A/i . 

Further examination of Fig. |51 indicates a notable ab- 
sence of any additional phase transition, as for example 
between isovector-isoscalar and possibly alpha-clustering 
phases. Clearly there are no phase separating lines far 
from the origin, i.e. at large At=o and At=i. Thus, a con- 
tinuous path between isovector and isoscalar phases does 
not involve a transitional behavior, unlike a transition 
from the normal state to any paired state. 

The point corresponding to the realistic nucleus. 
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T=0 
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1, is located in the transitional 
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region from the normal to isovector pairing phase on the 
side of the paired state. On the other hand, the isoscalar 
pairing matrix elements are to be increased by a factor of 
three to bring the system onto a border between the nor- 
mal phase and quasideuteron T = pairing coherence, 
this seen best from Fig. El Our result, that in the N = Z 
nucleus the T = 1 pair coherence is enhanced in contrast 
to T = 0, is supported by other findings 18, .25j. 

Let us summarize the results for our shell model ex- 
amples. In '*^Ca we are able to unambiguously identify 
the pairing phase transition. This nucleus would be in 
the paired state if no interactions of non-pairing type had 
been present. The incoherent residual interactions in re- 
ality put the ground state below the phase transition line 
although in the region of large fluctuations preceding the 
onset of well developed pairing. This explains earlier re- 
sults [i^: the failure of the BCS theory and large pair- 
ing correlation energy found in the exact solution. We 
have also found that the change in non-pairing interac- 
tions does not lead to a phase transition, however it influ- 
ences substantially the normal-to-pairing critical region. 
Weak residual incoherent interactions smear the phase 
transition as they facilitate the pair excitation. Strong 
non-pairing interactions create complexity in dynamics, 
and a strong phase transition occurs when such system 
is "cooled" and the superconducting order is established. 

Considering the ^''Mg nucleus we observed a clear 
phase transition between the normal and superconduct- 
ing state. We identified the regions where two different 
types of pairing, isovector and isoscalar, exist. The ^"^Mg 
nucleus is found, in agreement with earlier results, to be 



in the isovector phase and within the critical region of 
the normal to superconducting transition. The isoscalar 
pairing should have been three times stronger compared 
to the realistic shell model forces in order to lead to the 
corresponding condensate. Within limited shell model 
space we found no evidence of a phase transition between 
different types of pairing or alpha-clustering. 

To conclude, in this work we suggested a new theoreti- 
cal technique for tracking and quantitatively studying the 
phase transitions in nuclei and, more generally, in meso- 
scopic systems. Invariant correlational entropy serves not 
only as an indicator of the degree of complexity along the 
spectrum "T] but as a powerful tool for identifying the re- 
gions of rapid changes of structure in response to varia- 
tions in a control parameter, here it was a strength of the 
interaction. A mesoscopic system is sufficiently large to 
display statistical regularities but still sufhciently small 
to allow one to probe, theoretically and experimentally, 
individual wave functions. Such a system significantly 
alters the features of phase transitions, and opens a door 
for "secondary" interactions to leave a noticeable foot- 
print. The ICE clearly shows the sensitive spots of the 
parameter space. 
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